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In this work, fundamental results for carrier statistics in graphene 2-dimensional sheets and 
nanoscale ribbons are derived. Though the behavior of intrinsic carrier densities in 2d graphene 
sheets is found to differ drastically from traditional semiconductors, very narrow (sub- 10 nm) rib- 
bons are found to be similar to traditional narrow-gap semiconductors. The quantum capacitance, 
an important parameter in the electrostatic design of devices, is derived for both 2d graphene sheets 
and nanoribbons. 
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Graphene, a 2-dimensional (2D) honeycomb structure 
of carbon atoms, has generated intense interest recently 
[1, 2, 3, 4]. It has been now demonstrated that nar- 
row graphene nanoscale ribbons (GNRs) posses bandgaps 
that are tuned by the ribbon width [3]. These properties, 
along with the good transport properties of carriers (high 
mobility, high Fermi velocity) suggest that it is possible 
that graphene may be used in the near future in high 
speed electronic devices [6, 7]. In spite of rapid advances 
in the study of transport properties of graphene, basic 
tools of semiconductor device design such as temperature 
dependent carrier statistics and electrostatic properties 
such as quantum capacitance remain unexplored. This 
work investigates these properties for both 2D sheets, 
and GNRs in a comparative fashion, and analytical re- 
sults for these quantities are presented. 

The dispersion of mobile 7r-electrons in graphene in the 
1st Brillouin Zone (BZ) is given by 



E(k) = shv F \k\ 



(1) 



where s = +1 is the conduction band (CB), and s = — 1 
is the valence band (VB), h is the reduced Planck's con- 
stant, vp ~ 10 8 cm/s is the Fermi velocity of carriers 

in graphene, and |k| = 



i-2 



ky is the wavevector of 

carriers in the 2D (x — y) plane of the graphene sheet 
[4, 5]. The point |k| = 0, referred to as the 'Dirac point', 
is a convenient choice for the reference of energy, thus 
E(\k\ = 0) = eV. Each k— point is two-fold spin degen- 
erate (g s = 2), and there are two valleys in the 1st BZ 
(the K & K' valleys), g v = 2. 

In an undoped layer of graphene in thermal equilib- 
rium, there are mobile electrons in the CB and holes in 
the VB, similar to the intrinsic carriers in a pure bulk 
semiconductor. To find the 2D sheet density of such in- 
trinsic carriers in graphene, the linear density of states 
(DOS) 



Pgr(E) 
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(2) 



is used to write the 2D electron gas sheet density in 
graphene as 



dEp gr (E)f(E), 



(3) 



where f(E) is the Fermi-Dirac distribution function given 
by f(E) = (l + cxp[(E - Epj/kT})- 1 , k being the Boltz- 
mann constant, T the absolute temperature, and Ep is 
the Fermi-level. With the aid of the dimensionless vari- 
ables u = E/kT and r\ — Ep/kT (s is the index for 
conduction/ valence bands), the electron density may be 
rewritten as 



7r nvp 

and the hole density is symmetric, given by 
tt nvp 



(4) 



(5) 



where %(?/) = 1/T(j + 1) J °° du ■ u j /(l + e n ~'') is the 
Fermi-Dirac integral with j = 1. 

Under thermal equilibrium and under no external per- 
turbation (no applied bias, no optical illumination), the 
Fermi level is unique, and moreover, it is exactly at the 
Dirac point (Ep = eV). Then, the intrinsic carrier con- 
centration in 2D graphene is given by 



P 



TT kT_ 2 
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which is dependent on only one material parameter - the 
Fermi velocity. The point to note is that the intrinsic 
sheet density of electrons/holes does not depend on tem- 
perature exponentially; it has a T 2 dependence, due to 
the absence of a bandgap, and the linear energy disper- 
sion. At room temperature, the intrinsic electron and 
hole sheet densities evaluate to ~ 9 x 10 10 cm~ 2 . 

The situation changes for nanoscale ribbons cut from 
infinite graphene sheets. Consider a GNR of width W . 
Current experimental evidence suggests no clear depen- 
dence of the bandgap of GNRs on the chirality [3] . Re- 
gardless, the results derived here remain applicable for 
GNRs with bandgaps. We make the assumption that 
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FIG. 1: Intrinsic sheet carrier concentrations in graphene 
sheets and nanoribbons. Wide GNRs and 2d graphene have 
non-Arrhenius dependence on temperature, which becomes 
increasingly Arrhenius-like as the ribbon width decreases. 




FIG. 2: 2D carrier concentration in graphene and GNRs of 
different widths as a function of the location of the Fermi level. 
Though narrow GNRs exhibit large charge modulation due to 
the existence of a gap, they become similar to 2D graphene 
sheets when the Fermi level is deep inside the bands. 



the electron and hole quasimomenta are isotropic in the 
graphene plane. By aligning the x— axis along the longi- 
tudinal direction of the ribbon, the electron wavevector 
in the y— direction is quantized by hard-wall boundary 
conditions to be k y = nir/W (n = ±1,±2, ...), and the 
energy dispersion relation (Equation 1) for the n th sub- 
band becomes 



E(n, k x ) = shvf 



fc2 + (— ) s 



(7) 



indicating that the CB (n > 0) & VB (n < 0) split into 
a number of ID subbands, indexed by n. It is obvious 
that breaking the symmetry of a graphene sheet by cut- 
ting out a ribbon opens up a band gap. For the isotropic 
case assumed here, the bandgap for a GNR of width W 
is given by E g = £(+1,0) - £(-1,0) = 2nhv F /W, de- 
pendent only on the Fermi velocity and the width of the 
GNR. The DOS relation for the n th ID subband is then 
given by 



pGNn(n,E) 



E 



irHvp ^/E 2 - El 



®{E-E n ), (8) 



where 0(...) is the Heaviside unit step function, and 
E n = nirhvp/W = nE g /2. This directly leads to a to- 
tal DOS pgnr(E) = Y Jn p GNR {n,E). The expression 
for the total DOS is the same for the CB and VB, and 
exhibits van-Hove singularities at energies E n from the 
Dirac point. The electron density as a result is given by 



n 



AkT 
nhvf 



E 

n>0 



S(x n ,r)), 



where x n — E n /kT, r\ — Ep/kT, and 

u du 

b{x,r]) 



Vu 2 - x 2 1 + e u - v ' 



(9) 



(10) 



The intrinsic carrier concentration in GNRs is obtained 
when rj = Ep = 0, i.e., the Fermi level is at mid-gap; 
this leads to — (AkT/Tthvp) S(x n , 0). For narrow 
GNRs, E g 3> kT , then one can use the approximation 
S(x,0) « xKi(x) where Ki(...) is the Bessel function 
of 1st order, and the asymptotic approximation of the 
Bessel function Ki(x) ~ y // T:/2xexp (—a?) for large x to 
write the intrinsic carrier density of GNRs as 



(11) 



For bandgaps well in excess of the thermal energy, it 
suffices to retain only the first term in the summation to 
recover the familiar dependence 




/ irkT 



W\ E a 



(12) 



This relation has to be used with caution when ex- 
perimentally extracting bandgaps from the slope of 
Arrhenius-like plots; it is applicable only when the 
bandgap is well in excess of the thermal energy. The 
ID carrier concentration of GNRs may be converted to 
an effective 2D sheet density by writing ri2d = ^id /W for 
comparing their properties with graphene, as is done in 
Figure 1 . This figure shows that the intrinsic carrier con- 
centrations in GNRs differs significantly from graphene 
only if the ribbon widths are below ~0.1 pm, and in- 
dicates when Arrhenius dependence of intrinsic carrier 
concentrations is valid. 

The carrier sheet density in graphene can be changed 
by an electrostatic gate voltage, and the on-state sheet 
densities can approach, and exceed those in conventional 
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FIG. 3: Left: Quantum capacitance of 2D graphene and a 
5nm GNR compared with the parallel-plate capacitance of 1 
nm SiC>2 & HfG*2. Right: 2D carrier density in a graphene 
sheet as a function of gate voltage for different oxide thick- 



appear in the GNR density due to van-Hove singularities, 
as indicated by arrows. 

An important quantity in the design of nanoscale de- 
vices is the quantum capacitance [8]. Writing the total 
charge in a graphene sheet with a local channel electro- 
static potential V c h as Q = q(p—n) where q is the electron 
charge, and using the definition of quantum capacitance 
Cq = dQ/dVch, one obtains for 2D graphene 



2q 2 kT 
ir(hvp) 2 



In [2(1 + cosh 



qVchs 
kT ' 



which under the condition qV c h 3> kT reduces to 
2 2 qVch 



C c 



7T (flVp) 2 



q Pgr(qV ch ). 



(14) 



(15) 



If the electrostatic capacitance formed between a gate 
electrode and the graphene layer is given by C ox — 
tox/tox, then the electron density in the graphene layer 
can be written as a function of the gate voltage as 



field-effect transistors (FETs) . If the Fermi level in a 2D 
graphene sheet is driven from the Dirac point to Ep = 
r]kT electrostatically by means of a gate voltage, then the 
electron density is given by n — 7^^1(77) /^(O), and the 
hole density by p = rj)/^si(0), leading to a mass- 

action law np — nf 3i(ry)9i(— 77)/3^(0). Similarly, if the 
local electrostatic potential in a GNR is tuned by a gate 
voltage such that the Fermi level is at Ep — r\kT, then 
the electron density is given by equation 9. For r\ > x 3> 
1, one can make the approximation (1 + exp {u — ~ 
0(77 — u) to rewrite the ID electron concentration as 



E 



EIQ(E F -E n ), (13) 



On the other hand, for a non-degenerate condition when 
the Fermi level is located inside the GNR bandgap, us- 
ing the approximation S(x,rf) ~ yJiTxj2 exp{r] — x) the 
electron concentration may be written as n ~ rite 11 , and 



similarly, for holes, p 



;e v , which is similar to tra- 



ditional semiconductors. Figure 2 shows the calculated 
exact 2D carrier concentrations in graphene and GNRs of 
different widths as a function of the location of the Fermi 
level (qV c h — Ep) at room temperature. Though narrow 
GNRs exhibit large charge modulation due to the exis- 
tence of a gap, they become similar to 2D graphene sheets 
when the Fermi level is deep inside the bands. Ripples 



" =tig- n Q { ,/l + 2— - 1), 
1 nn 



(16) 



where uq = C ox VG/q is the traditional carrier density 
one would obtain by neglecting the quantum capacitance, 
and no — \(C ox h,vpjq 2 ) 2 , which arises solely due to the 
quantum capacitance. Applying the same method to find 
the quantum capacitance (per unit width) of GNRs, one 
obtains for the condition r) > x ^> 1, 



Cc 



Aq 2 



E 



= 0(ry - x n ) = q 2 pGNR{ri). 



(17) 

The quantum capacitance of 2D graphene and GNR is 
plotted in Figure 3 (left), and compared with the oxide 
gate capacitance of 1 nm Si02 & Hf0 2 . Figure 3 (right) 
shows the carrier density dependence in 2D graphene on 
the gate voltage (Eq. 16) for different Si02 thicknesses. 
Gate modulation of the charge is strong but non-linear 
for very thin t ox since Cq ~ C ox under that condition. 
The field-effect becomes weak but increasingly linear as 
t ox is increased since C ox « Cq under that condition. 

The results presented here would prove useful for the 
design of electronic devices using graphene sheets and 
GNRs. The authors would like to thank G. Snider for 
useful discussions, and one of the authors (D.J.) thanks 
an NSF CAREER award for financial support. 
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